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A gauge-invariant quantum theory of the Friedmann-Robertson- Walker (FRW) universe with dust 
is studied in terms of the Ashtekar variables. We use the reduced phase space quantization which has 
following advantages: (i) fundamental variables are all gauge invariant, (ii) there exists a physical 
time evolution of gauge-invariant quantities, so that the problem of time is absent and (iii) the 
reduced phase space can be quantized in the same manner as in ordinary quantum mechanics. In 
the FRW model, the dynamical components of the Ashtekar variables are given by a single quantity 
p and its conjugate momentum, where p is related to the scale factor a as a oc i/]pj and its sign gives 
the orientation of triads. We solve a scattering problem in terms of ingoing and outgoing energy 
eigenstates. We show that the incident wave is reflected in rate 1/4 and transmitted in rate 3/4 at 
the classical singularity p — 0. Analyzing the dynamics of a wave packet, we show that the classical 
initial singularity is replaced by a big bounce in quantum theory. A possible interpretation of the 
result is that the wave function of the universe has been in a superposition of states representing 
right-handed and left-handed systems before the big bounce. 



PACS numbers: 04.60.-m, 04.60.Ds, 98.80.Qc, 04.20.Fy 



I. INTRODUCTION 



The initial singularity is inevitable in most of the classi- 
cal cosmological models. One of the motivations of quan- 
tum cosmology is to shed light on quantum nature of 
the singularity. However, there exists potential problems 
that have not been completely resolved yet. A problem 
is about what should be interpreted as observables in 
classical and quantum gravity 0, HJ . A canonical formu- 
lation of general relativity (GR) is a constrained system 
with hrst-class constraints in which the spacetime diffco- 
morphisms are interpreted as gauge transformations. In 
gauge theories, only gauge-invariant quantities are ob- 
servables. However, there are technical and conceptual 
difficulties in the realization of the idea especially in GR 
(see, e.g., [3| and references therein). In many works, 
gauge-variant quantities are used as observables. This 
issue must be seriously considered especially in quantum 
gravity because it is substantially related to the problem 
of time [1] . 

Recently, a gauge-invariant quantum theory of the 
Friedmann-Robertson- Walker (FRW) universe with dust 
has been proposed |5] . There a variable associated with 
the dust serves as an internal time @. The so-called re- 
lational formalism 0, II] has been used to construct the 
classical reduced phase space spanned by gauge-invariant 
quantities, and then the system has been quantized in the 
same manner as in ordinary quantum mechanics. In the 
method, the phase space is reduced at the classical level 
and then quantized, so that the procedure is completely 
different from the Dirac quantization [9| where one quan- 
tizes the unconstrained phase space at first and then im- 
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poses the quantum constraints on physical states. Al- 
though the full gauge-invariant quantum theory of grav- 
ity has not been obtained at present, the method gives a 
possible resolution to the problem of time and observable 
at least in the context of minisuperspace models. 

It has also been shown in [5] that the initial singu- 
larity is replaced by a big bounce. In the analyses, the 
scale factor a itself has been chosen as one of the canon- 
ical variables and its range has been restricted to a > 
to avoid the negative scale factor. One may think that 
the big bounce is caused by the boundary conditions on 
wave functions due to the restriction, that is, the result 
depends on the variable one choose. In fact, since canon- 
ical formulations that were equivalent classically yield 
different quantum theories, it is quite interesting to in- 
vestigate the presence of the big bounce by using other 
well-motivated variables with unrestricted range. 

In this paper, we shall construct and analyze the 
gauge- invariant quantum theory of the flat FRW universe 
with the dust in terms of the Ashtekar variables [T3, E[ ■ 
The Ashtekar formulation has been frequently used in 
quantum gravity (e.g. [12]) and, in particular, is em- 
ployed as a startin g p oint of the quantization in loop 
quantum gravity [si. Il3j|. In the Ashtekar formulation, a 
SU(2) connection A l a and a densitized triad Ef constitute 
a canonical pair where, in the FRW model, the dynami- 
cal component of the triad is given by a single quantity 
p such that a oc y/\p\. The range of p is the whole real 
line, which is doubled from that of a, and the sign of 
p determines an orientation of the triad. We use the 
reduced phase space quantization method as in the pre- 
vious work, which is different from loop quantization. In 
the present system, since the Hamiltonian becomes sin- 
gular at p = 0, we discuss carefully the self-adjointness of 
the Hamiltonian and the space of wave functions. In view 
of the quantization, the Ashtekar variables are shown to 
be quite natural because there are no need to impose any 
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boundary conditions on wave functions at the singularity. 

We shall then analyze the obtained quantum theory 
to reveal the quantum nature of the initial singularity. 
First, we solve a scattering problem in terms of ingoing 
and outgoing energy eigenstates. The analysis shows that 
a part of the incident wave is reflected and the rest is 
transmitted at p = 0. Second, we consider the motion 
of a wave packet and evaluate the expectation value of 
the scale factor. It is shown that the expectation value 
has a non-zero minimum, that is, the initial singularity 
is replaced by a big bounce in quantum theory. The 
remarkable point is that the big bounce mixes the states 
representing right-handed and left-handed systems. One 
of the obtained scenarios is that if the present universe 
is in the right-handed state, the past universe before the 
big bounce was in a superposition of states representing 
right-handed and left-handed systems. 

The organization of the present paper is as follows. In 
Sec. UH we introduce the classical reduced phase space of 
the FRW universe with dust. In Sec. lIIIl we quantize the 
reduced system and obtain the Schrodinger equation. In 
Sec. II VI we analyze the obtained quantum theory to see 
the dynamics of the universe. In Sec. [Vj we conclude the 
paper, and Sec. |VI]is devoted to discussion. 

In this paper we adopt the unit in which c = 1. 



II. CLASSICAL SYSTEM 

In this section, we shall construct the reduced phase 
space of the classical FRW cosmology. 

First, we present the classical system of the FRW 
universe with dust in terms of the Ashtekar vari- 
ables. Second, we briefly review the relational formal- 
ism, emphasizing that the relational formalism for a de- 
parametrized system enables one to explicitly construct 
gauge-invariant quantities and a physical Hamiltonian 
which generates the time evolution thereof. Third, we 
construct the reduced phase space of the FRW universe 
with dust by using the relational formalism. 



A. Classical phase space 

In the Ashtekar formulation, the variables (A l a ,Ef) 
form a canonically conjugate pair where A l a is a SU(2) 
connection and Ef is an orthonormal triad with density 
weight 1. We here consider the flat FRW universe which 
is described by the metric 



-dt 2 + a 2 {t){dx 2 + dy 2 +dz 2 ). 



(1) 



In the flat FRW model, the Ashtekar variables can be 
written in terms of only one independent components c 
and p, 



4 = E?=p(t)X?, 



(2) 



where u/ are bases of left invariant one-forms and X, 
are invariant vector fields dual to the one-forms. These 
variables have relations to the scale factor a such that 



1 

c = sgn(p)— d, 



(3) 



where 7 is the so-called Barbero-Immirzi parameter [Til . 
Il4j |. N is the lapse function and the dot denotes the 
derivative with respect to t. Note that, while the scale 
factor is restricted to be nonnegative, p ranges over the 
entire real line, carrying an orientation of triads deter- 
mined by the sign of p. 

Since the three-space integral, defined as V := J dx 3 , 
diverges in the flat FRW model with R 3 topology, one 
must somehow get rid of this divergence. There are at 
least two ways to treat this issue. The first is to con- 
sider a compact universe. The second is to introduce a 
finite cell and restrict the range of the integral to this 
cell. The physical interpretation of the resulting theory 
becomes clearer in the first case because one can interpret 
Ipp/ 2 ) defined below as the length of the shortest non- 
trivial loops of the universe. Thus, we take the first op- 
tion and in particular we only consider the case of three- 
dimensional torus, where we take a cube of coordinate 
range < x,y,z < Vz , and identify the opposite faces. 

If we define new variables as p := V^p and c := V~5c, 
the gravitational action is written as 



--"grav 



dt 



— pc + NH glav 
«7 



(4) 



where k = 8nG and the Hamiltonian constraint H. 
takes the form 



grav 



K7 Z 



(5) 



In terms of the new variables, the Poisson brackets are 
independent of V: 



K.7 



(6) 



Turning to matter couplings, we consider the dust in- 
troduced by Brown and Kuchaf Q . The dust is dynami- 
cally coupled to gravity and the action of the dust in the 
FRW spacetime takes the form 



S, 



dust 



dt 



P T T - NH, 



dust 



(7) 



Here, T is the proper time measured along the particle 
flow lines when the equations of motion hold, Pt is its 
conjugate momentum and 



H 



dust 



= Pt 



(8) 



In particular, it can be shown by solving the equations 
of motion that T coincides with the cosmological proper 
time in isotropic and homogeneous models. It serves as a 
natural time variable in the construction of the reduced 
phase space. 
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Combining these results, we now have the total action 
for gravity plus the dust, 



St„t = / dt 



—pc + P T T - NH tot 
K7 



(9) 



where 



Htot — flgrav + -ffdust 



j(?V\p\ + p T = 0- (10) 



The form of the constraint (fT0| is called deparametrized 
form, which has a special feature as will be seen in the 
discussion below. 



gauge orbit 




B. Relational formalism 



FIG. 1: The action of the map Ci H on phase space points and 
the gauge- invariance of Op(y). 



We here give a brief review of the relational formalism, 
which is a method to construct gauge-invariant quantities 
in constrained systems (see @, 3 for details). From now 
on, we assume that the system has only one constraint 
for simplicity. 

The key observation of the relational formalism to de- 
fine gauge-invariant quantities is as follows. Take two 
gauge-variant functions F and T on the phase space, and 
choose one of the functions T as a clock. Then, the value 
of F at T = t is gauge-invariant even if F and T them- 
selves are gauge variant. That is, a selected function T 
serves as a clock and the relation between T and other 
variables is interpreted as time evolution. 

Let us now move to the mathematical definition of 
the gauge-invariant quantities in the relational formal- 
ism. Suppose a phase space has a 2n-dimension (n > 2), 
and there are canonical coordinates (q a ,p a , a — 1, - ■ ■ ,n) 
such that {q a ,Pb} = Sb- We- will denote a first-class con- 
straint by H and a phase space point by y — (q a ,p a ). 
Under the gauge transformation generated by H, a point 
y is mapped to y i-> a t H (y), where t is a gauge parameter. 
That is, a^iy) is a gauge flow generated by H starting 
from y. 

Then we can define the gauge- invariant quantity T F (y) 

as 



Op(y) ■= F(a t H (y))\ T(a t i(y)> 



(11) 



Here, F(a t H (y)) can be written as a series F(a t H (y)) = 

YZU%{ H > F h»)(v)' where { H ^ F }(0) ■■= F and 
{H,F} (n+1) := {H,{H,F} {n) }. The definition dTTJ) 
gives a manifestly gauge-invariant quantity because 
Op(y) is constant on each gauge orbit. Indeed, as illus- 
trated in FigQ] if y and y' are on the same gauge orbit, 
there is t' such that oic(y) — a c(y') an( l 



O^y) = F (aUy))\ T(aUy)hT 
= F(a t r(y')) 



OW)- (12) 



A constraint equation H = is said to be of de- 
parametrized form if it is written as 



H(q a ,T,p a ,P T ) =P + h(q a ,p a ) = 



(13) 



with some phase space coordinates {q a , T;p a , Pt}- 
In the deparametrized theories, the reduced phase 
space is spanned by the gauge-invariant quantities 
(Oqa (y), Op a (y)) associated with q a and p a with the sim- 
ple symplectic structure 



{o r g «(y),o; b (y)} = si 



(14) 



The physical Hamiltonian H is obtained by replacing q a 
and p a in h(q a ,p a ) with T qa {y) and OJ a (y), 

H (OJa (y), T Pa (y)) := h (OJ, (y), T Pa (y)) . (15) 

The Hamiltonian generates the time evolution of the 
gauge-invariant quantity associated with a function F 
which depends only on q a and p a : 



dO T F {y) 
dr 



= {H,OUv)}- 



(16) 



As we have seen, the relational formalism enables one 
to construct gauge-invariant quantities and extract their 
time evolution in the deparametrized case. Therefore, the 
method provides us a possible resolution to the problem 
of time and observable. 



C. Reduced phase space of 
Friedmann-Robertson- Walker universe with dust 

Let us now construct the reduced phase space of the 
FRW universe with dust. As mentioned above, if the 
constraint is written in the deparametrized form, one 
can construct the reduced phase space coordinatized by 
gauge-invariant quantities and extract their time evolu- 
tion by using the relational formalism. 
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In the present case, it is natural to choose the function 
T as the clock variable. Then, the reduced phase space is 
coordinatized by the gauge-invariant quantities C(r) := 
Ol(y) and P(t) := Op(y) associated with c and p with 
very simple symplectic structure 



{C(t),P(t)} = ^. 



(17) 



Moreover, we can obtain the physical Hamiltonian H p h_ ys 
by replacing c and p in _ff grav (c, p) with C and P, 



phys 



K7 Z 



rC(r) 2 #(i. 



(18) 



The Hamiltonian (|18|) generates time evolution of a gauge 
invariant quantity F as 



dF 



{H phys ,F}; F = F(C(t),P(t)). (19) 



So far, we have constructed the reduced phase space 
of the FRW universe with dust. We will consider the 
quantization of the system in the next section. 



III. QUANTIZATION 

In this section, we shall quantize the system on the 
reduced phase space obtained in the previous section. 

The quantization procedure is straightforward as in the 
ordinary quantum mechanics because we need not to deal 
with the constraint. Now the physical variables are oper- 
ators and the Poisson bracket {•, •} is replaced with the 
commutation relation (l/ih)[», •]. Thus (fTT)) becomes the 
canonical commutation relation 



[C,P} = 



(20) 



The evolution of the state \& of the universe is described 
by the Schrodinger equation 



ih 



a7 



^phys ^ i 



(21) 



where H phys := H phys (C,P). 

Although the prescription above formally gives a quan- 
tum theory of the universe, the theory is not well de- 
fined or uniquely determined until we specify the oper- 
ator ordering, the Hilbert space and the the domain of 
the operator i/ph ys which makes H p h ys self-adjoint. Self- 
adjointness, beyond Hermiticity, of H p h ys is substantial 
in quantum mechanics because it is necessary for unique 
development of the quantum state and the conservation 
of probability. 

Let us choose the ordinary Schrodinger representation 
in which the operators P and C , respectively, act on a 
wave function ^(P) in the following way: 

P*(P) = P*(P), (7* ( P) = ^p^ffi. (22) 



As a concrete example, we choose the following opera- 
tor ordering for the Hamiltonian, 



^ntivs 



(23) 



Then the Schrodinger equation (|2"Tj) takes the simplest 
form, 



iH 



d 2 ^ 

Qp-2- 



(24) 



Since the present Hamiltonian is different from the or- 
dinary kinematical term, we choose the Hilbert space as 
H = L 2 (R, \P\~zdP) the set of square integrable func- 
tions of P with respect to the measure \P\~?dP in order 
to make the Hamiltonian (|23|) Hermitian up to surface 
term. We shall rewrite the Schrodinger equation (|24"|) in 

2 

dimensionless variables x := Pj 3 lp\ and f := r/tpi, 



df 



dx 2 



(25) 



where lp\ and ipi are the Planck length and the Planck 
time, respectively. 

Though the Hamiltonian -ff p h ys is somewhat singular 
at the origin x = 0, it is indeed self-adjoint in %, which 
we shall show in Appendix |XJ Note that we do not need 
any special boundary conditions for the wave function ^ 
at x — 0. 



IV. DYNAMICS OF THE UNIVERSE 

In this section, we shall investigate the obtained one- 
dimensional quantum system to analyze the quantum na- 
ture of the initial singularity. First, we solve a scattering 
problem in terms of ingoing and outgoing energy eigen- 
states. Second, we analyze the motion of a wave packet 
and evaluate the expectation value of the scale factor. 
Third, we construct a special analytic wave function and 
investigate its time evolution. 



A. Reflection and transmission rates 

We here calculate the reflection and transmission rates 
at the origin x — 0. From the Schrodinger equation (|25l) . 
the energy eigenvalue equation reads 



H pbys ij:(x) = y/\x\ 



d 2 iP{x) 
dx 2 



Etp{x). 



(26) 



We here note that w\x\d /dx is a negative definite self- 



adjoint operator on L 2 (l 



'dx). The negative defi- 



niteness, (ip\Hp^ ya \tp) < for all ip, can be easily shown 
by integration by parts, and it is consistent with the fact 
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that the gravitational energy never takes positive val- 
ues in classical theory. Then, we can rewrite the energy 
eigenvalue equation ([26)) as 



dV(s) 
dx 2 



(27) 



To obtain the eigenf unctions, we make use of the fact 
that the equation (|27[) can be transformed into Bessel's 
differential equation. Let us first find the solution of 
Eq. (|2"T|) for x > 0. In terms of <fi — ip/y/x and 
]E\x§, Eq. (J2U yields 



dz 2 

The solution is written as 



deb , , 4 



)0 = O. 



(z) = A 1 J 2/3 (z) + A 2 J_ 2/3 (z) 



(2), 



'2/3 



2/3^ 



(28) 

(29) 
(30) 



where Aj., A 2 , Bi and B 2 are arbitrary constants, J± 2 /3 
are the Bessel functions and Hy' 3 are the Hankel func- 
tions. Thus, the eigenfunctions for x > take the form 



Ax J. 



2/3 



A 2 J_ 



2/3 



BxH. 



(i) 

2/3 



E\xi 



B 2 H, 



(2) 



2/3 



E\x* 
(31) 

£|a; 
(32) 



For later use, we here explicitly give the definition of 

(-1)" (z\ 2n + v 



the special functions; J v {z) := X^^Lo ' 



(1,2) 



(z) 



and ff, 



function N u (z) is defined as N v {z) 



n\T{n+v+l) 

J v [z) ± iN v (z) where the Neumann 

J v (z) cos(z^7r) — J— v (z) 

sin(i>7r) 

To obtain the energy eigenfunctions in the whole x, we 
connect the functions for x > and x < at the origin 
with the conditions ip(—0) = tp{+0) and dip/dx(—0) = 
dip/dx(+0). We have 



ijj(x) 



A 1 sgn(x)J 2/3 



E\\x\* 



-A 2 J^ 



2/3 



Ettx]* 



(33) 



where the first term is an odd function and the second 
is an even function. Here we should note that these en- 
ergy eigenvalue functions do not lie in the Hilbert space 
L 2 (R, dx). This situation is similar to the system of 
a free particle in quantum mechanics and one can expand 
a physical state in terms of the energy eigenstates. 

Let us now move to the calculation of the reflection 
and transmission rates. We recall here that the Han- 

(1 21 

kel functions have the asymptotic forms Hi (z) ~ 
_2_ e ±«[z— ("+5)f ] for large z. Thus they are the ana- 



logues of the exponential functions e 



t IX 



kx, in the 



B„ 



x=0 



FIG. 2: The conception of the incident, reflected and trans- 
mitted waves. 



free particle case. We consider the following stationary 
state ip(x); (i) the incident wave propagates in the posi- 
tive x direction and its amplitude is set to one, and (ii) 
there are no waves that propagate in the negative x di- 
rection for x — y oo, that is, 



i>{x) =B Iei ^/\x~\H^l (ly/\E\\x\l 



x\H, 



(2) 

2/3 ^3 

r(i) 



E\\x\ 4 for x < 0, 



(34) 



■»»-vix\H?/ 3 [^V\E\\x\*) forx>0, (35) 



where B Te [ and B trans are the amplitudes of the reflected 
and transmitted waves, respectively (See Fig. [2|). We 
can calculate the amplitudes by imposing the conditions 
-0(-O) = V>(+0) and dtp/dx(-Q) = dip/dx(+0). These 
conditions yield 



B Te 



z-V3' 



BtlD 



V3 

-V3" 



(36) 



The difference between the phases of the reflected and 
the transmitted waves is — ir/2 because B tlans / B rc [ = 
V3e~ l f . 

In the present system, the conservation law for proba- 
bility density -M?p(x, f ) + -§^j(x, f ) = holds where the 
energy density p(x, f ) and its flow j(x, f ) are defined as 



p(x,f) 



\^{x,f)\ 2 



(37) 



j(x,f) := i 



dx 



Ox 



(38) 



Then, one can define the reflection rate R and the trans- 
mission rate T as 



R = 



3r{x) 



T = 



3t{x) 



•7/0*0 



(39) 



where jj, jn and jj>, respectively, are the flow of the 
probability density for the incident wave, the reflected 
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wave and the transmitted wave. We can obtain the values 
of R and T at the origin as follows: 



R\x=o = 
TL=o = 



3r(0) 



3i (P) 
Jt(0) 



= \B Tei \ 2 = 



— -Btrs 



1 

4' 



(40) 
(41) 



Note that J2| x= o and T| x =o do not depend on the energy 
eigenvalue E. So far, it has been shown that the incident 
wave reflects in rate 1/4 and transmits in rate 3/4. That 
is, although we do not impose any boundary conditions 
at the origin, neither a complete reflection nor a complete 
transmission can occur in the present model. 




T [tj - 006 



P [(3/&cH£] 



FIG. 3: The absolute value of the wave function is plotted as 
a function of r and P. 



B. Time evolution of wave packet and scale factor 

Let us now analyze the dynamics of a wave packet. 
The procedure is as follows. First, we prepare an initial 
wave packet ^{x, 0) at some nonzero x. Then, we nu- 
merically evolve it backward in time by the Schrodinger 
equation (|25p and evaluate the expectation value of |x| 
as a function of the internal time f. Here we consider |x| 
because both the positive and negative x correspond to 
the universe of the same size with different orientation of 
triads. 

The numerical methods used here are the forth-order 
Runge-Kutta method in the time integration and the 
midpoint difference method for the spatial differentiation. 
However, we need a special care for several points near 
the origin because of the singular nature of the Hamil- 
tonian. The midpoint difference method does not give 
a good approximation at the origin, because there the 
right hand side of Eq. (|25p incorrectly becomes zero all 
the time. Thus, we use the least square method to fit the 
wave function near the origin. We consider the eigenfunc- 
tion expansion of the wave function and take the first five 
terms in the series as the fitting function f(x): 

f(x) =ci+ c 2 x + c 3 \x\? + c 4 sgn(x)\x\^ + c 5 \x\ 3 . (42) 

Then the right hand side of Eq. (f^5|) can be evaluated as 

r, — 7<P f 3 15 „ . .3 

VH^ = -c 3 + —c 4 x + 6c 5 \x\i . (43) 

For simplicity, we here choose the initial wave function 
as a Gaussian wave packet 

0) = Co exp (- (3: ~ o ,2 0)2 - ikox^J , (44) 

where Co is the normalization constant. Although the 
other choices are possible, one can obtain qualitatively 
same results irrespective of the choice. In the numerical 
calculation below, we set the initial values as io = 5, 
fc = 20 and a — 0.5. Fig. [3] shows the absolute value 
of the wave function as a function of P and t, and the 
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FIG. 4: The expectation values of P is plotted as a functions 
of the time r. 



expectation value of \P\ is plotted as a function of the 
time t in Fig 2] 

We can see from Fig. [3] that, as is expected from the 
previous discussion, a part of the wave packet is reflected 
and the rest part is transmitted at the origin. We here 
remind that the sign of x determines an orientation of tri- 
ads, which correspond to a right-handed and left-handed 
systems respectively. Thus, the result indicates that if 
the present state of the universe is in a right-handed sys- 
tem, the past state is in superposition of the states of a 
right-handed and left-handed systems. 

This calculation allows another interpretation because 
^*(x, — f ) is the time-reversal solution of the Schrodinger 
equation. If we take complex conjugate of the initial state 
(|4"4")l and set f = — f' , we can consider the initial wave 
packet as a state in the past and the above calculation as 
forward evolution in time. That is, if a past state of the 
universe has been in a right-handed system, the present 
state is in superposition of the states of a right-handed 
and left-handed systems. 

As for the expectation value of \P\, Fig. 2] indicates 
that the expectation value never goes to zero and bounces 
at a nonzero minimum. That is, the initial singularity is 
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replaced by a big bounce in the present model. The result 
is consistent with the previous work [5J, where one does 
not need any boundary conditions at the origin. The 
difference is that the wave function in this model is a 
superposition of the states representing the right-handed 
and left-handed systems. 



C. Special analytic wave function 

Let us construct a special analytic wave function. We 
superpose the odd functions in the energy eigenstate ([33| : 



i, i\x) = sgn(x)^\x~\J 2/3 (jy/\E\\x\) 



(45) 



The superposition of the states is written in the form 



whe 




dd. 



C(e)e^ £ %J 2/3 (|x|le)d6, (46) 



E\ and C is a superposition coefficient. 



_ 8 £ 2/3 e -ae 2 



We choose the function C(e) to be C(e) — g , 
where a is an arbitrary positive constant. Then we can 

use the formula J °° e~ ax2 x v+1 J u (bx)dx = ( 2 a)^-+ i e ~^° 
where Re(a) > and Re(v) > -1 (e.g., [IH). The inte- 
gration yields the wave function 



\&(x,f) 



(2/3(r))S 



-exp 



4/8(f) 



(47) 



where /3(f) := a — i^-f. The wave function (|4T|) has 
double peaks, which represents a equal- weight superpo- 
sition of the right-handed and left-handed states. 

The expectation value of \x\ also can be calculated as 



(N) 



(*IMI*> 



r(f) 



16. 



(48) 



The solution again shows that the initial singularity is 
avoided. For large f , the scale factor a is evaluated as 



KT3. 



(49) 



Thus, they are in good agreement with the classical tra- 
jectories when the universe is sufficiently large. 



V. CONCLUSIONS 

A gauge- invariant quantum theory of the flat FRW uni- 
verse with dust has been studied in terms of the Ashtekar 
variables. In the FRW model, the dynamical component 
of the Ashtekar variables are given by a single quantity p 



and its conjugate momentum c, where p is related to the 
scale factor a as a oc \/\p\, and the sign of p determines 
an orientation of the triad. We have first constructed the 
classical reduced phase space of the system by using the 
relational formalism and then have quantized the reduced 
system. The advantages of the quantization method are 
as follows: (i) fundamental variables P and C are gauge- 
invariant quantities associated with the Ashtekar vari- 
ables p and c, (ii) a natural time evolution of the gauge- 
invariant quantities exists, so that the problem of time 
and observable is resolved at least in the deparametrized 
theory and (iii) the reduced phase space can be quantized 
in the same manner as in ordinary quantum mechanics 
because there are no constraints in the reduced phase 
space. In the quantization, we have rigorously discussed 
and have shown the self-adjointness of the Hamiltonian 

H phys = --^j\/\P~\C 2 which is singular at P = 0. There 
we need not to impose any boundary conditions on wave 
functions at the singularity because the range of P is not 
restricted to be finite as opposed to the scale factor. 

In the obtained quantum theory, we have analyzed 
the quantum nature of the initial singularity. We have 
first considered a scattering problem by using the Han- 
kel functions. We have shown that the incident wave is 
reflected in rate 1/4 and transmitted in rate 3/4 at the 
origin P = which was the initial singularity in clas- 
sical theory. Then, we have investigated the motion of 
a wave function of the universe. In the present pa- 
per, we have considered different two scenarios in which 
the initial wave functions ^(P, 0) have been chosen re- 
spectively as a wave packet peaked at some positive P 
and an anti-symmetric wave which has double peaks. In 
both cases, it is shown that the expectation value of P 
has a non-zero minimum, that is, the initial singularity 
is replaced by a big bounce in quantum theory. The 
interpretation for the case of a wave packet is twofold 
depending on backward or forward evolution. From the 
backward evolution, it follows that if the present state of 
the universe is in a right-handed system, the past state 
has been in a superposition of the states of a right-handed 
and left-handed systems, and from the forward evolution, 
it follows that if a past state of the universe has been in 
a right-handed system, the present state is in a superpo- 
sition of the states of right-handed and left-handed sys- 
tems. In the case of the anti-symmetric wave function, 
the universe has remained in a equal-weight superposi- 
tion of the right-handed and left-handed states all the 
time. 



VI. DISCUSSION 

The solution scattered at the origin P — may be 
related to the existing imbalance of the parity of the uni- 
verse. 

A possibility is that the universe after a big bounce 
is a superposition of the right-handed and left-handed 
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states, and one of them is chosen by destruction of the 
superposition. If one takes into account the degrees of 
freedom of local gravitational field and other fields, the 
superposition may be considered as a quantum cat state. 
The state is vulnerable to gobservation" , or other small 
perturbations, which breaks symmetry of parity. 

Another possibility is that the scattering solutions pre- 
sented in this paper correspond to local phenomena in the 
universe, namely, collapses and formation of black holes. 
For example, if one considers situation of a Oppenheimer- 
Snyder collapse where the inner spacetime is described by 
a FRW spacetime, one expects that inside the black hole 
there is a "big bounce" beyond which there is a universe 
of the opposite parity. Then there may be regions in 
the universe where there are elementary particles with 
opposite parity from those that we presently know. 



continuity (any absolute continuous function is an inte- 
gral of some L 1 function). 

Symmetricity of i? p hy S is simple. One can easily 
show by partial integration that -ffphys is Hermitian, i.e., 
(<S>,H phys V) = (H phys <f>,$), for G D(H phys ). Then 
-ffphys is symmetric because Z3(i? p h ys ) is dense in H. 

Let us find the expression and the domain of H^ hys . 

Let A be the restriction of H phys onD(A) := G H\~$> G 
Cq°(1R)}, where Cq°(R) is the set of smooth functions on 
R with compact support. For f ef(A) and Seflf^), 
wc have 



f (a1$) = (At$, ) = (</>, Ail>) = j dx<$><S>". 

J ■y \x\ J 

(A2) 
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Appendix A: Self-adjointness of -Hphys 



Eq. (|A2[) implies that the weak second derivative Z? 2 $ of 
<f> is given by D 2 $ = A^/y/~\x\. bmce At$ g H and / G 



Let us show that H p \ iys = \/\x\-^z is self-adjoint in 
the Hilbert space U = L 2 (R, \x\-2dx) where the inner 
product is given by ($,5') := J H dxQ^ / y/\x\. Difference 
from the case of standard quantum mechanical Hamilto- 
nian H = ~-^ + V{x) on L 2 (R) with well-behaved V(x) 
(e.g. [lj|) is that we need some care on the point x = 0. 
However, it turns out to cause no harm. 

We define the domain of H p h ys by 

D{H phys ) = {* G H\f G C\R),^' is locally 

absolutely continuous, iJphys^ G "H}, 

(Al) 

where a prime denotes the derivative, and the second 
differentiation in H p h ys is in the sense of local absolute 



U implies f/\x\^ 4 G L 2 (R,dx), we have L> 2 $ = g/\x\^ 4 
with g G L 2 (R,dx). 

We can show that D 2 ^ is locally L a (R,dx) with 1 < 
a < 4/3. By Holder's inequality, we have, for norms 
on any compact set K, Wg/lx] 1 / 4 ^ = ||(<7/|£| 1//4 ) a || 1 ^ 
H- Q/4 IL - ILall"lllx|-^/ 4 ||} /9 where q = 



1, lllxl-^/ 4 !!! hence 



115 112/q IIFI ' ll« = H5II2 II 
2/(2 -a). Since qa/4 < 1, |||x| -~< 
||5/k| 1/4 || Q is finite - Thus is locally L a (R,dx). 

Therefore $ is locally W 2 - a (R,dx). 

By Sobolev's embedding theorem, this implies $ is in 
c i,2-i/a( R y In p ar ticular, $ is in C^R), i.e., $ is 
differentiable everywhere and has a continuous deriva- 
tive This $' must be locally W 1 ' a (R, dx) because $ 
was locally W 2,a (R, dx). Since a > 1, is also locally 
VF 1,:L (R, dx), namely, $' is locally absolutely continuous. 
Thus any <3? G D(A) is differentiable and has a locally 
absolutely continuous 

It follows that any $ G D(H^ hys ) has a locally ab- 
solutely continuous derivative for we have D(A^) D 
D ( H t hys \ from D(A) C D(H phys ). Therefore ffj hys ap- 

plies as H^ hys ^ = y^x^", where the second differentia- 
tion is in the sense of local absolute continuity. We also 
have D(Hp hys ) — Z?(i7 p hy S ) because A is Hermitian on 

D (Hphys). 

As a result, H p h ys is self- adjoint on its domain (|Alj) . 
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